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We consider a discrete model that describes a locally regulated 
spatial population with mortality selection. This model was studied 
in parallel by Bolker and Pacala and Dieckmann, Law and Murrell. 

We first generalize this model by adding spatial dependence. Then we 
give a pathwise description in terms of Poisson point measures. We 
show that different normalizations may lead to different macroscopic 
approximations of this model. The first approximation is determin¬ 
istic and gives a rigorous sense to the number density. The second 
approximation is a superprocess previously studied by Etheridge. Fi¬ 
nally, we study in specific cases the long time behavior of the system 
and of its deterministic approximation. 


1. Introduction. We consider a spatial ecological system that consists of 
motionless individuals (such as plants). Individuals are characterized by their 
location. We assume that each plant produces seeds at a given rate. When a 
seed is born, it immediately disperses from its mother and becomes a mature 
plant. We also assume that plants are subjected to mortality selection. That 
is, each plant dies at a rate that depends on the local population density. 
All these events occur randomly in continuous time. This model was intro¬ 
duced by Bolker and Pacala [2] and Dieckmann and Law [9] . To study the 
system, Bolker and Pacala derived approximations for the time evolution of 
the moments (mean and spatial covariance) of the population distribution. 
In the present article, we wish to give a rigorous definition of the underlying 
microscopic stochastic process and rewrite rigorously the moment equations 
of [2], then to derive some tractable macroscopic approximations, and finally 
to study the long time behavior of the stochastic process and its approxi- 
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mations. Unfortunately, we obtained only partial results concerning the last 
point. 

In Section 2, we describe the Bolker-Pacala-Dieckmann-Law (BPDL) 
process in detail. In fact, we generalize the model slightly by adding a spa¬ 
tial dependence in all the rates. Then we give a pathwise representation 
of the system in terms of Poisson point measures. We also produce a nu¬ 
merical algorithm to simulate the BPDL process. Section 3 is devoted to 
existence and uniqueness. We also show some martingale properties of the 
BPDL process. In Section 4, we find the mean equation that Bolker and 
Pacala [2] intuitively obtained. We also give a rigorous sense to the covari¬ 
ance terms formally defined in [2] or [9], [4] and [10]. Section 5 is concerned 
with macroscopic approximations of the BPDL process. We first show that, 
conveniently normalized, the BPDL process converges to the solution of a 
deterministic nonlinear integrodifferential equation. We propose this as a 
rigorous interpretation of the density number, often introduced by biologists 
without a proper definition. We also show that with another normalization, 
the BPDL process converges to the superprocess version of the BPDL model 
introduced and studied by Etheridge [6]. We give partial results about ex¬ 
tinction and survival for the BPDL process in Section 6. In Section 7, we 
study the convergence to equilibrium of the deterministic approximation. 
We obtain only some partial results. We next show that in the detailed bal¬ 
ance case to be specihed later on, there exists a nontrival steady state for 
the BPDL process. We conclude the article with some simulations. 

2. The model. Let us first describe the model in detail. 

2.1. Definition of the parameters and heuristics. The plants are sup¬ 
posed to be motionless and characterized by their spatial location. We as¬ 
sume that the spatial domain is the closure X of an open connected subset 
X of for some d>l. We denote by Mp{X) [resp. V{X)] the set of finite 
nonnegative measures (resp. probability measures) on X. Let also Xi be the 
subset of Mp{X) that consists of all finite point measures: 



( 2 . 1 ) 


Here and below, 5x denotes the Dirac mass at x. For any m = ^xi ^ Xi, 
any measurable function / on X, we set {m, f) = f dm = f{xi). 

Notation 2.1. For all x in .T, we introduce the following quantities: 

(i) fi{x) G [0,oo) is the rate of “intrinsic” death of plants located at x, 

(ii) 7 (x) G [0,oo) is the rate of seed production of plants located at x, 
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(iii) D{x,dz) is the dispersion law of the seeds of plants located at x. It 
is assumed to satisfy, for each x G X, 

/ _D{x,dz) = l and / 

JzeM.‘^,x+z£A: Jzi 

(iv) a{x) G [0, oo) is the rate of interaction of plants located at x. 

(v) For X, y in X, U{x,y) = U{y,x) G [0,oo) is the competition kernel. 


_ D{x, dz) = 0. 


The competition kernel U{x,y) describes the strength of competition be¬ 
tween plants located at x and y. 

We aim to study the stochastic process ut, taking its values in Ai and 
describing the distribution of plants at time t. We write 

i(t) 

(2.2) ft = dxi, 

i=l 

where I{t) gN stands for the number of plants alive at time t and ,..., 
describe their locations (in X). The supposed dynamics for this population 
can be roughly summarized as follows: 


(i) At time f = 0, we have a (possibly random) distribution uq G Ai. 

(ii) Each plant (located at some x G X) has three independent expo¬ 
nential clocks: a seed production clock with parameter 7 (x), a natural death 
clock with parameter fi{x) and a competition mortality clock with parameter 

a{x)El^lu{x,Xl). 

(iii) If one of the two death clocks of a plant rings, then this plant disap¬ 
pears. 

(iv) If the seed production clock of a plant (located at some x G X) rings, 
then it produces a seed. This seed immediately becomes a mature plant. Its 
location is given by y = x + z, where z is randomly chosen according to the 
dispersion law D{x,dz). 


In [2], 7 , /i, a and D were assumed to be space-independent. Our general¬ 
ization might allow us to take into account external effects such as landscape, 
resource distribution and so forth. Note also that assuming that all these 
clocks are exponentially distributed allows us to reset all the clocks to 0 each 
time one clock rings. 

We wish to describe the system by the evolution in time of the empirical 
measure Uf. More precisely, we are looking for an Al-valued Markov process 
(^t)t>o with infinitesimal generator L, defined for a large class of functions 
(p from Ai into K, for all v G Ai, by 


(2.3) 


L(l){u)= u{dx) D{x,dz)[(f){u + 5x+z) - 
Jx JRrf 

+ J^u{dx)[(j){i' - 4 ) - 4>{i^)]^^y{x) + a{x) J^i^{dy)U{x,y) 
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The first term is linear (in v) and describes the seed production and dispersal 
phenomenon. The second term is nonlinear and describes death due to age 
or competition. This infinitesimal generator can be compared with formula 
(3) in [2], page 182. 

2.2. Description in terms of Poisson measures. We now give a pathwise 
description of the Al-valued stochastic process To this end, we use 

Poisson point measures. For the sake of simplicity, we assume that the spatial 
dependence of all the parameters is bounded in some sense. 

Assumption A. There exist some constants a, 7 and p, such that, for 
all X G A, 

(2.4) a{x)<a, 7 (x)< 7 , Ai(x) < fi. 

There exist a constant C > 0 and a probability density D on such 
that, for all X G A, 

(2.5) D{x,dz) = D{x, z) dz with. D{x, z) < CD{z). 

The competition kernel U is bounded by some constant U. 

We also introduce the following notation. 

Notation 2 . 2 . Let N* = N\{0}. Let H = {H\ ..., ...) :M ^ (M'^)^* 

be defined by 

(2.6) = (x^(i),...,X^(„),0,...,0,... ), 

where Xo-(i) ^ ^ for some arbitrary order ^ on (one may, e.g., 

choose the lexicographic order). 

This function H allows us to overcome the following (purely notational) 
problem: Assume that a population of plants is described by a point measure 
u G . Choosing a plant uniformly among all plants consists of choosing i 
uniformly in { 1 ,..., (u, 1 )}, and then choosing the plant number i (from the 
arbitrary order point of view). The location of such a plant is thus Lf*(u). 

Notation 2.3. We consider the path space T C B([0, 00 ), Mi?(A)) de¬ 
fined by 

1 /Vf > 0,ut G Af, and 30 = fo < < ^2 < • • •, 

^ lim„f„ = 00 and ut = ut. Vf G [fj,fi+i) 

Note that for {r't)t>o ^ T, and f > 0 we can define in the following way; 
If t ^ = ut; while if t = ti for some z > 1, vt- = • 
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We now introduce the probabilistic objects we need. 

Definition 2.4. Let {Q,T,P) be a (sufficiently large) probability space. 
On this space, we consider the following four independent random elements: 

(i) an Af-valued random variable uq (the initial distribution); 

(ii) a Poisson point measure N{ds,di,dz,d6) on [0,oo) x N* x x [0,1], 
with intensity measure yds iJ2k>i ^k{di)){CD{z) dz) dO (the seed production 
Poisson measure); 

(hi) a Poisson point measure M{ds,di,d6) on [ 0 , 00 ) x N* x [0,1], with 
intensity measure fids {J2k>idk{di)) dO (the “intrinsic” death Poisson mea¬ 
sure); 

(iv) a Poisson point measure Q{ds,di,dj,d6,d9') on [0,oo) x N* x N* x 
[0,1] X [0,1], with intensity measure Ua ds {J2k>i dk{di)){J2k>i dk{dj)) dO dO' 
(the “competition” mortality Poisson measure). 

We also consider the canonical filtration {Pt)t>o generated by these pro¬ 
cesses. 

We finally write the BPDL model in terms of these stochastic objects. 

Definition 2.5. Admit Assumption A. A (d^i)i>o-adapted stochastic 
process u = {vt)t>o that belongs a.s. to T will be called a BPDL process if 
a.s., for all t>0, 

LLi i(i<(-.-.i))V(..-)+.) 

X N{ds,di,dz,d6) 

(2-8) X*i) ^ii<('^^-d)}^Hfius-)^{e<{i^{Hfius-)))/mMids,di,d9) 

~i LLi I 

X ^{e<{a{Hfius-)))/{a)}Qids, di, dj, d9, d9'). 

Although the formula looks complicated, the principle is very simple. The 
indicator functions that involve 9 and 9' are related to the rates and appear 
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when the parameters depend on the space variable x. In the case where the 
rates are constant (studied in [2]), all the integrals and indicator functions 
that involve 6 may be cancelled. 

Let us now show that if u solves (2.8), then it follows the dynamics in 
which we are interested. 

Proposition 2.6. Admit Assumption A. Consider a solution 
to (2.8). Then {i't)t>o is a Markov process. Its infinitesimal generator L is 
defined for all bounded and measurable maps 4>:A4^'R, all v € A4, by (2.3). 
In particular, the law of {nt)t>o does not depend on the chosen order {see 
Notation 2.2). 

Proof. The fact that {vt)t>o is a Markov process is classical. Let us 
now consider a function (f as in the statement. Recall that with our nota¬ 
tion, = Recall also that L(/)(r'o) = dtE[(f){vt)]t=o- A simple 

computation, using the fact that a.s. 4>{vt) = 4>{t^o) + J2s<t[4>i^s- + — 

Ps-}) — 4>{i^s-)], shows that 





X N{ds, di, dz, dO) 










X Q{ds,di,dj,d9,d6'). 


Taking expectations, we obtain 



h iCD{z) 
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+ / dsE 

Jo 


+ / dsE 






{Us,l) (Us,l) 

E'“E E 

i=l j=l 


[Us- - dHi{u,_)) - (t){l^s-)] 


U{H\vs-),W{vs-))a{H\us-)) 


U 


a 


= E[(j){vQ)\-\- [ dsE [_i's{dx) [ dz'y{x)D{x, z) 

Jo Ux JRd- 

X + < 5 (x+ 2 )) - 4>{^s)] 


+ / ds E 


I i^s{dx)[4>{i^s - 4 ) - (Pi’^s)] 

Ux 

X ^d{x) + aix) j_ Vs{dy)U{x, y) 
Differentiating this expression at t = 0 leads to (2.3). □ 


2.3. About simulation. This pathwise definition of the BPDL process 
leads to the following simulation algorithm: 


Step 0. Simulate the initial state tzq and set Tq = 0. 

Step 1. Compute the total eventuate, given by m(0) = mi(0)+ 7712 ( 0 ) + 
m 3 ( 0 ), with 

(2.9) mi(0) = C'7(po,l), "^ 2 ( 0 ) =/7(7/o, 1), msiO) = aUiizoAf■ 

Simulate Si exponentially distributed, with parameter m(0), and set Ti = 
To + S*!. Set vt = t'o for all t <Ti. Choose whether to go to Step 1.1, 1.2 or 
1.3 with probability mi(0)/m(0), m2(0)/m(0) and m3(0)/m(0). 

Step 1.1. Choose i uniformly in {1,..., {uq, 1 )}. Choose 2 : G according 
to the law D{z) dz. With probability 1 — {uo), z))/{^CD{z)), 

do nothing (i.e., set vti = 1 ^ 0 ); else, add a new plant at the location f7®(po) + 2 
(i.e., set UTi = no + ^(m{uo)+z))- 

Step 1.2. Choose i uniformly in {1,..., (z/q, 1)}. With probability 1 — 
{fi{W{nQ)))/fl, do nothing (i.e., set = z^o); else, remove the zth plant 
(i.e., set z/Ti = no - 
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Step 1.3. Choose i and j uniformly in {1, ..., {uq, 1)}^. With probability 
1 — {vQ))a{W{vQ)))/Ua, do nothing (i.e., set vti = else, 

remove the ith plant (i.e., set vti =^o — ^h^{uo))- 


Step 2. Compute the total event rate, given by m(Ti) = mi(Ti) + 
m 2 {Ti) +m 3 {Ti), with 

mi{Ti) =C'7(i/Ti,1), 


( 2 . 10 ) 


"12(^1) =/r(pTi,l), 

msiTi) = aU {nTiAf ■ 


Simulate S 2 exponentially distributed, with parameter m(Ti), and set T 2 = 
Ti + 82 - Set vt = for all t£ [Ti,T 2 [ and so forth. 

3. Existence and first properties. We now show existence, uniqueness 
and some moment estimates for the BPDL process. 

Theorem 3.1. (i) Admit Assumption A and that E{{i2q,1)) < 00 . Then 

there exists a unique BPDL process the sense of Definition 2.5. 

(ii) If furthermore, for some p>l, E{{vq, 1)^) < 00 , then for any T < 00 , 



(3.1) 


El sup 

\t£[0,T] 


Proof. We hrst prove (ii). Consider thus a BPDL process We 

introduce for each n the stopping time = inf{t > 0, {vt, 1) > n}. Then a 
simple computation using Assumption A shows that, neglecting the nonpos¬ 
itive death terms, 

sup {vs,lY 


sg[0,tAT„] 


< 





X N{ds, di, dz, dO) 



X N{ds, di, dz, dO) 


(3.2) 
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for some constant Cp. Taking expectations, we thus obtain, the value of Cp 
changing from line to line: 


E 


sup 

:fo,tAT„l / 


- s€[0,tAT, 

(3.3) <Cp + CpE 
< Cp + CpE 


rt/\Tn 


ds'jC dzD{z)[{iys_,l) + Ws-Ay] 


ds[l + {l^sATr,Cf] 


The Gronwall lemma allows us to conclude that for any T < oo, there exists 
a constant Cp^T, not dependent on n, such that 


(3.4) 


E 


sup <Cp,T- 

O.TATr7l / 


tG[0,TAT, 


First, we deduce that tends a.s. to infinity. Indeed, if not, we can find a 
To < oo such that sto = P^sup^Tn < Tq) > 0. This would imply that for all 
n, T(supjg[o,ToArn](^u 1)^) ^ ^TafP-, which contradicts (3.4). We may let n 
go to infinity in (3.4) thanks to the Fatou lemma. This leads to (3.1). 

Point (i) is a consequence of point (ii). Indeed, we can, for example, build 
the solution (t't)t>o using the simulation algorithm previously described, 
and choosing the rates and acceptance-rejection according to the Poisson 
measures N, M and Q. We have to check only that the sequence of (effective 
or fictitious) jump instants T„ goes a.s. to inhnity as n tends to infinity, and 
this follows from (3.1) with p=l. Uniqueness also holds, since we have no 
choice in the construction. □ 

We now prove that if there is at most one plant at each location at time 
t = 0, then this also holds for all t>0. 

Proposition 3.2. Assume Assumption A and that E{{vq, 1)) < oo. As¬ 
sume also that a.s., sup,j,g_:p ^ 1- Consider the Bolker-Pacala process 

(^'t)t>o- Then for all t > 0, a.s., 


(3.5) / nt{dx)nt{{x}) = {nt,l) 

Jx 


that is, sup t't({x}) < 1. 

x&X 


Proof. Consider the nonnegative function cj) dehned on Ai by = 
A p((ix)p({x}) — (p, 1). Then note that a.s. <?i>(po) = 0 and that for any v G 
Ai, any x G suppz/, (j){n — 5x) — < 0- Consider, for each n > 1, the 

stopping time = inf{t > 0, {nt, 1) > n}. A simple computation allows us to 
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obtain, for all t>0, all n > 1, 


(3.6) 


<0 + E 



D{x, dz)'y{x) 




We easily check, using that u is atomic, that the right-hand side term identi¬ 
cally vanishes, since D{x, dz) has a density. Hence, a.s., 4>{i^tATn) = 0- Thanks 
to (3.1) with p = 1, Tn a.s. grows to infinity with n, which concludes the 
proof. □ 


We carry on with a property that concerns the absolute continuity of 
the expectation of ut- For v a random measure, we define the deterministic 
measure E{u) by {E{u), f) = E{{u, f)). 


Proposition 3.3. Accept Assumtion A, that E[{uq,1)] < oo and that 
E^uq) admits a density no with respect to the Lebesgue measure. Consider 
the BPDL process Then for all t > 0, E{vt) has a density np, for all 

measurable nonnegative functions f on X, E[{i't,f)] = X? dx f{x)nt{x). 


Proof. Consider a Borel set A of with Lebesgue measure zero. Con¬ 
sider also, for each n > 1, the stopping time = inf{t > 0, {ut, 1) > ra}. A 
simple computation allows ns to obtain, for all t>0, all n > 1, 


T’[(PtAr„,lA)] =£^((p0,1a)) 


(3.7) 


+ E 


- E 


( ds [ n.idxMx) [ 

\Jo JX JK 


ds _iAsidx)''y{x) / dz D{x, z)1a{x + z) 


ptATn 


ds _i's{dx)lA{x) 
Jx 


X (h^{x) + a{x) J_Us{dy)U{x,y) 


By assumption, the first term on the right-hand side is zero. The second term 
is also zero, since for any x e A, f^d dz 1a{x + z)D{x, z) = 0. The third term 
is of course nonpositive. Hence for each n, E{{vtATn-, 1a)) is nonpositive and 
thus zero. Thanks to (3.1) with p = 1, Tn a.s. grows to infinity with n, which 
concludes the proof. □ 


We hnally give some martingale properties of the process 
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Proposition 3.4. Admit Assumption A and that for some p > 2, EKuq, 1)^] < 
oo. Consider the BPDL process o,nd recall that L is defined by (2.3). 

(i) For all measurable functions 4> from M into M such that for some 
constant C, for all v G A4, |(/>(p)| + \L(l){v)\ < C(1 + (z/, 1)1*), the process 


(3.8) 


fiut) - - [ dsL4>{us) 

Jo 


is a cadlag -{J-t)t>o-f^o,'<'l'l'ngale starting from 0. 

(ii) Point (i) applies to any measurable (p satisfying \4>{v)\ < C'(1 + (p, 

(iii) Point (i) applies to any function 4>{h') = fu^fY, with 0 < q < p — 1 
and with f bounded and measurable on X. 

(iv) For any bounded and measurable function f on X, the process 


(3.9) 


Xl/= (utj) - {uqJ) - [ ds [ Us{dx)j{x) [ dzD{x,z)f{x + z) 

Jo Jx JRd 


+ [ ds [_iys{dx)f{x) 

Jo Jx 


/x(x) + a(x) _i^sidy)U{x,y) 


lx 


is a cadlag -martingale starting from 0 with {predictable) quadratic vari¬ 
ation 


(3.10) 


{Mf)t= [ ds [ iys{dx)\j{x) [ dz f{x + z)D{x,z) 
Jo Jx I JK'* 


+ f{x) 


p.{x) + a{x) _Us{dy)U{x,y) 


lx 


Proof. First of all, note that point (i) is immediate thanks to Propo¬ 
sition 2.6 and (3.1). Points (ii) and (iii) follow from a straightforward com¬ 
putation using (2.3). To prove (iv), we first assume that Pl[(z/o,l)^] < oo. 
We apply (i) with ^{ly) = {iy,f)- This yields that is a martingale. To 
compute its bracket, we hrst apply (i) with (p{iy) = {iy,fY obtain that 


(3.11) 


{i^tjf - {lyojf 

— ds _iys{dx)j{x) / dzD{x,z) 

Jo Jx JR<i 

X [f{x + z) + ‘^f{x + z){i^s,f)] 

- [ ds [ iyfidx){f{x)-2f{x){iys,f)} 

Jo Jx 


ia{x) + a{x) i'fidy)U{x,y) 

Jx 
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is a martingale. Then we apply the Ito formula to compute (z/t,/)^ from 
(3.9). We deduce that 

Wtjf - Wojf 


(3.12) 


- [ ds [ i's{dx)'Y{x) [ dzD{x,z)2f{x + z){us,f) 
Jo Jx JmJ- 

rt 


+ [ ds [ Vs{dx)2f{x){h's,f) 
Jo Jx 


fi{x) + a{x) Vs{dy)U{x,y) 

Jx 


- {Mf)t 


is a martingale. Comparing (3.11) and (3.12) leads to (3.10). The extension 
to the case where only ill[(t'o, 1)^] < oo is straightforward since, even in this 
case, < oo thanks to (3.1) with p = 2. □ 


4. On the the BPDL moment eqnations. We now wish to give a sense 
to the mean moment equation given in [2], formula (6). Note that in the 
biology literature, one may be confused by the notation between the discrete 
measure ut, its expectation £'(t't) [defined by {E{iyt),f) = E{{i't, f))] and a 
measure with density nt{x) of which the definition is not clear. Following [2] 
in this section we use the next assumption. 


Assumption B. The spatial domain is A = M'^. All the parameters a, 
7 , // and D of the model are independent of x. Moreover, the (bounded) 
competition kernel U{x,y) has the form U(x — y), and both dispersal and 
competition kernels are symmetric probability distribution functions, that 
is, D{z) = D{—z), U{x — y) = U{y — x) and /j^d dzD{z) = fj^d dzU{z) = 1. 


We moreover assume that 1)^) < oo and that there is at most one 

plant at each location at time t = 0. So (3.1) with p = 1 holds and we can 
define, for each time t € [0,T], 

(4.1) n(t) = 


Using Proposition 3.4(iv) with / = 1 and taking expectations in (3.9), we 
obtain 


(4.2) 


E((i^i,l))=E((iyo,l)) + 



(j-y)E((iy^,l)) 


— a 




iZs{dx)us{dy)U {x 
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Hence, 


(4.3) 


n(t) = n( 0 ) + (7 —/i) [ dsn{s) — a [ dsE( f Us{dx)U{0)us{{x}) 
Jo Jo \Jm.‘^ 


— a dsE 
Jo 


Us{dx)i^s{dy)l{xj^y}U{x - y) . 


However, thanks to Proposition 3.2, we know that for all s > 0, /j^d Vsidx)!! (0) x 
i's{{x}) = U{0){i'sA) ■ We thus obtain 

n(t) = n(0) + (7 —/r — a?7(0)) [ dsn{s) 

(4.4) ° 

— a I dsE 


[ dsE( [ Us{dx)iys{dy)ls^^yxU{x-y) 
Jo \JR<^xR‘i 


Let us now explain the covariance term used by Bolker and Pacala. Writing 


aE 


(4.5) 


^ Us{dx)us{dy)l{^^y}U{x - y) 


= aE[ I ns{dx){us{dy)-n{s)dy)l{xj^y}U{x-y)]+an'^{s), 
we obtain, from (4.4), 

n(t) = n( 0 ) + (7 —/i — al7(0)) [ dsn{s) — a [ dsn^{s) 

Jo Jo 


(4.6) 


— a ds E 
Jo 


us{dx){us{dy) - n(s) dy)li^^yjU(x - y) 


Following the terminology of Bolker and Pacala, we define a covariance mea¬ 
sure Ct on for each time t. Let T-y denote the translation by the vector 
—y. We set 

(4.7) Ct{dr) = e(^J ^{r^o}t^t o rZy{dr) 0 ntidy)^ - n^(t) dr. 


In other words, the covariance measure is defined for each measurable bounded 
function cj) with compact support in by 


Ct{dr)(j){r) 


(4.8) 


= E 

= E 




vt{dx)vt{dy)l^xi^y}(l){x - y^ -n^{t) dr4i{r) 


_ vt{dx){vt{dy) - n{t) dy)l{^^y}(t){x - y) 
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By using this notation, we obtain the mean equation obtained by Bolker and 
Pacala ([2], formula (6), page 183), with a rigorous sense for the quadratic 
term: 

(4.9) — n(t )(7 — /i — an{t)) — aU{0)n{t) — a [ Ct{dr)U{r). 

at jRd 

Let us finally remark that we are also able to derive an evolution equation for 
the covariance measure. In other words, we can write differential equations 
solved by Ct{dr)(p{r) for all measurable bounded functions (j) on (we, 
however, do not obtain the same equation as in [2]). Of course moments of 
higher order are involved in such equations. So a remaining issue is to find 
reasonable moment closures as developed in [4]. These closures are, at the 
moment, not rigorously justified. 

5. Infinite particle approximations. Our aim in this section is to describe 
the effect of two different normalizations on the BPDL process. In both cases, 
we make the initial number of plants grow to infinity. We first consider 
the case where the birth and death rates are unchanged. We show that 
the random measure {nt)t>o tends to a deterministic measure {^t)t>o and 
solution of a nonlinear integrodifferential equation. 

In addition, the second normalization consists of accelerating the rates in 
a convenient way. Then {nt)t>o converges to a superprocess (W)t>o- This 
measure-valued process was introduced by Etheridge [6], who called it the 
superprocess version of the Bolker-Paeala model. 

Let us first consider the most general situation. 

Notation 5.1. For each n G N*, we consider a set of parameters {finj'Jn, Oin,Un, Dn) 
as in Notation 2.1, that satisfy for each n. Assumption A and consider an 
initial condition Vq £ Ai. Then, we denote by {nf’)t>o the BPDL process 
(see Dehnition 2.5) with the corresponding coefficients. We consider the 
subset Al” of Mf{X) defined by 

(5.1) = 

We finally consider the cadlag Af"-valued Markov process {Xf')t>o defined 
hyX^ = lvf. 

The generator of (A")t>o is then given, for any measurable map (j) from 
A4" into R, by 

L'^(f){v)=n [ v{dx) [ dzjn{x)Dn{x,z) (l)(iy + -dx+z 

, , Jx JRd. _ \ n 

(5.2) 
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+ n J_u{dx)^^Hnix) + nanix) j_ v{dy)Un{x, y) 


Mv - 5x ] - (t){v) 


n 


Indeed, the generator L” of (r'")t>o is given by (2.3), replacing (/r, 7 , a, U, D) 
by {iin,1n,an,Un,Dn). Hence, 

(5.3) LXz.) = dtE,[cj){X^)]t=o = dtEnu[cl^{i^?/n)]t=o = L^rini^), 

where (/)"'(^) = The conclusion follows from a straightforward com¬ 

putation. We now restate Proposition 3.4 for the renormalized model. 


Lemma 5 . 2 . Let n>l be fixed and consider the process {X^)t>o defined 
in Notation 5 . 1 . Assume that for some p>2, E[{Xq, 1)^] < 00 . 

(i) For all measurable functions from At"" into M such that for some 
constant C, for all v G A4", \(t>{n) \ -|- \ L'^(j){v)\ < C(1 -|- {v, 1)^), the process 

( 5 . 4 ) fiiXf) - fiiXS) - f dsL^fiiX^) 

Jo 

is a cadlag -martingale starting from 0. 

(ii) Point (i) applies to any measurable satisfying \4>{n)\ < C{l + {u, 

(hi) Point (i) applies to any function 4>{n) = {i>,f)^, with 0 < q < p — 1 

and with f bounded and measurable on JH. 

(iv) For any f bounded and measurable on X, the process 

Mfi’^ = {XfJ)-{XS,f) 

( 5 . 5 ) - [ ds f_Xg{dx) [ dz'yn{x)Dn{x,z)f{x + z) 

Jo Jx Jmd 

+ Jo ds J^X^{dx)^^y.nix) nan{x) J^X^{dy)Un{x,y)^ f{x) 

is a cadlag -martingale with {predictable) quadratic variation 

= - / ds [ X'^{dx) [ dz'^n{x)Dn{x,z)f{x + z) 
n Jo Jx JRd 

( 5 . 6 ) + - [^ ds f X^{dx) 

n Jo Jx 

X Jj,n{x)+ nan{x) j_X'f{dy)Un{x,y)'^f{x). 

We endow Mp{X) with the weak topology. 
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5.1. Convergence to a nonlinear integrodifferential equation. Let us now 
consider the mean-field approximating case in which the initial number of 
particles n tends to infinity, and the parameters of seed production and in¬ 
trinsic death stay unchanged, whereas the mortality competition parameter 
tends to zero as We show that the BPDL process can be approximated 
by a deterministic nonlinear integrodifferential equation. This might be a 
better deterministic way to describe the model than the moment equations 
of [2]. In particular, it allows us to deal with space-dependent parameters. 

Assumption Cl. 

1. The initial conditions converge in law and for the weak topology on 

Mp (A) to some deterministic hnite measure Mp{X), and sup„ E ( {Xq , 1)^) < -|-oo. 

2. There exist some continuous nonnegative functions 0,7 and /r on A, 
bounded by 0,7 and / 2 , such that ^n{x) = 7 (x), fin{x) = ti{x) and an(x) = 
a{x)/n. 

3. There exists a bounded nonnegative symmetric continuous function U on 
A X X bounded by U such that Un{x,y) = U{x,y). 

4. There exists a continuous nonnegative function D on X that satisfies, 
for each x £ X, x+zex dzD{x, z) = 1, D{x, z) = 0 as soon as x + z ^ 

X and such that D[x, z) < CD{z) for a constant C > 0 and a probability 
density D on M'^. We set Dn{x^z) = D{x,z). 

The hrst assertion of Assumption Cl is satished, for example, if Aq = 
n ^ 1=1 where the random variables Z® are independent, with law ^o- In 
this case, the number n can be seen as the volume of particles at initial time, 
and the limit of Xf — n^t xnay give a rigorous sense to the number density. 

Theorem 5.3. Admit Assumption Cl, and consider the sequence of pro¬ 
cesses A®® defined in Notation 5.1. Then for all T > 0, the sequence (A”) 
converges in law, in D([0, T], Mj?(A)), to a deterministic continuous function 
ift)t>o £ C{[0,T],Mp{X)). This measure-valued function is the unique 
solution, satisfying sup^gjQ 2 ’] (6; 1) < 00 , of the integrodifferential equation 
written in its weak form: for all bounded and measurable functions f from 
X into M, 

(6,/) = (?o,/) + / ds f is{dx)y{x) [ dzD{x,z)f{x +z) 

(5.7) _ 

- ds J_f,s{dx)f{x)^^fi{x) +a{x) J_f,s{dy)U{x,y)y 

Note that the link between (2.8) and (5.7) is the same as the link between 
the continuous-time binary Galton-Watson process with birth rate 7 and 
death rate y,, and the deterministic differential equation f'{t) = (7 — y)f{t). 
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Proof. We divide the proof into several steps. Let us fix T > 0. 


Step 3. Let us first show the uniqueness for equation (5.7). We consider 
two solutions (6)t>o and (6)t>o of (5.7) that satisfy supjg[o^r] (6 + 6,1) = 
At < + 00 . We consider the variation norm dehned for /ii and ^2 in Mf{X) 
by 


(5.8) 11/21-^211= sup K/Xi -/i 2 ,/)|. 

/eL°°(A),||/||oo<i 


Then we consider some bounded and measurable function / dehned on X 
such that II/II 00 < 1 and we obtain 


1 ( 6 - 6 ,/)! 


< ds 
Jo 


lx 


(5.9) 


+ / ds 
Jo 

+ f ds 
Jo 


[is{dx) -£,s{dx)\ 

X ( 7 (x) / dzD{x,z)f{x + z)-i^{x)f{x) 

V Jr<^ 

_[^s{dx) - ^s{dx)]a{x)f{x) f_^s{dy)U{x,y) 


L [^s{dy) - ^s{dy)] is{dx)a{x)f{x)U (x, y) 


lx 


lx 


However, since ||/||oo < 1 for all x G X, 


while 


and 


7 (x) / dzD{x,z)f{x + z)-y.{x)f{x) 
jRd- 


<7 + h, 


oi{x)f{x) / is{dy)U{x,y) 

Jx 


< odJ At 


[ is{dx)a{x)f{x)U{x,y) 
Jx 


< all At- 


We deduce that 

(5.10) 1(6-6,/)l < [7 + h + 2dt/^T] / ds||6-6ll- 

Jo 

Taking the supremum over all functions / such that ||/||oo 6 1 and using 
the Gronwall lemma, we hnally deduce that for all t <T, ||^i — 611 = 0- 
Uniqueness holds. 
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Step 4. Let us prove some moment estimates. By Assumption Cl, it is 
easy to check that, for all T > 0, 

(5.11) sup Li ( sup (Ai”,!)^ j <+ 00 . 

Vtg[0,T] / 

Indeed, recalling that we can prove, following line by line the 

proof of Theorem 3. l(ii) withp = 3, that Li[supig[o,r] (^” 1 1)^] — 1)^], 

noting that the constant Ct does not depend on n. We easily conclude using 
part 1 of Assumption Cl. 


Step 5. We first endow Mf{X) with the vague topology, the extension 
to the weak topology being handled in Step 6. To show the tightness of 
the sequence of laws Q"'= C{X^) in V{O{[0,T], Mf{X))), it suffices, fol¬ 
lowing [15], to show that for any continuous bounded function / on X, the 
sequence of laws of the processes (A”,/) is tight in D([0,T],M). To this end, 
we use the Aldous criterion [1] and the Rebolledo criterion (see [7]). We have 
to show 


(5.12) 


supL;[ sup \{X^J)\) <QO, 
n VtG[o,r] / 


and the tightness, respectively, of the laws of the martingale part and of 
the drift part of the semimartingales (A”,/). Since / is bounded, (5.12) is 
a consequence of (5.11). Let us thus consider a couple {S,S') of stopping 
times satisfying a.s. 0 < S < S' < S + 5 <T. Using Lemma 5.2, we get 


U[|M"/-M^’^|] 

(5.13) < < E[{M'''-^)s+5 - (M^’-^)5]^/2 


< E 


rS+S 


{j + p + m) I ds((A,^l) + (A^,l)2) 


1/2 


<cVs, 


where the last inequality comes from (5.11). The finite variation part of 
{Xs'J) - {XsJ) is bounded by 


rS+S 


(5.14) 


ds[{-f + p){X^,l) + aU{X^,lf] 

<601+ sup (Af,l)2 
se[o,r] 


Hence, formula (5.11) allows us to conclude that the sequence Q'' = £(A”) 
is tight. 
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Step 6. Let us now denote by Q the limiting law of a subsequence of 
Q”. We still denote this subsequence by Q”. Let X = {Xt)t>Q a process with 
law Q. We remark that by construction, almost surely, 

(5.15) sup sup |(Xf,/) - (Xf-,/)| < 1/n. 

This implies that the process X is a.s. strongly continuous. 


Step 7. Let us now check that a.s. the process X is the unique solution 
of (5.7). Thanks to (5.11), it satisfies supig[o,T](^u 1) < +oo a.s. for each 
T. Standard density arguments show that it suffices to check that X solves 
(5.7) for all / G Cb{X) and all t > 0. Let thus / G Cb{X) and t > 0 be hxed. 
For u G C'([0, oo),Mp{X)), denote 

^t(i^) = {^tj) - WoJ) 

(5.16) — [ ds f_i'a{dx)^{x) 

Jo Jx 

+ j ds l_Va{dx)f{x){^i.{x)+a{x) fvsidy)U{x,y) 

J 0 J PC L J PC 

We have to show that 


dz D{x, z)f{x + z) 


(5.17) ^Q[|Ti(X)|]=0. 

However, Lemma 5.2 and Assumption Cl imply that for each n, 

(5.18) M”’7 = Ti(A’^). 


A straightforward computation using Lemma 5.2, Assumption Cl and (5.11) 
shows that 


(5.19) E[\M^ 


nJi2] _ 


c 


= E[{M^’J\a < 

n 


f ds{l + {X^,If} 
Jo 


< 


C 




n 


which goes to 0 as n tends to infinity. On the other hand, since X is a.s. 
strongly continuous, since / is continuous and thanks to Assumption Cl, the 
function is a.s. continuous at X. Furthermore, for any v G ©([O, T],Mf{X)), 


(5.20) |^t(p)| < sup (1 + (z^s,!)^). 

sG[0,t] 


Hence using (5.11), we see that the sequence {^t{X'^))n is uniformly inte- 
grable and thus 

(5.21) limE(|Tt(A”)|) = .F(|Tt(A)|). 

Associating (5.18), (5.19) and (5.21), we conclude that (5.17) holds. 
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Step 8. The previous steps imply that the sequence (X”) converges to 
^ in ]D){[0,T],Mp{X)), where Mf{X) is endowed with the vague topology. 
To extend the result to the case where Mf{X) is endowed with the weak 
topology, we use a criterion proved in [12]: Since the limiting process is 
continuous, it suffices to prove that the sequence ((X"',l)) converges to 
(^,1) in law, in B([0,T],T). We may of course apply Step 5 with / = 1, 
which concludes the proof. □ 

Proposition 5.4. Assume that in Mf{X) has a density with respect 
to the Lebesgue measure. Consider the associated solution {^t)t>o io (5.7). 
Then for every t>0, the finite measure has a density with respect to the 
Lebesgue measure. 

Proof. The proof is similar to that of Proposition 3.3. We consider a 
Borel subset ^4 of T with measure zero. We apply (5.7) with / = 1^. The 
right-hand side expression is equal to 0 since the first term is zero by hypoth¬ 
esis, the second one is zero since for all x, fg^d dz lx+zeAD{x, z) = 0, and the 
last term is nonpositive. 

□ 

Remark 5.5. (i) Equation (5.7) is the weak form of, for all x G T, t > 0, 


(5.22) 


dtf.t{x)= I dyf,t{y)l{y)D{y,x-y) 

Jx 

- y{x)it{x) - a{x)ft{x) _ dyf,t{y)U{x,y). 

Jx 


(ii) Assume now that X = that the competition kernel is of the form 
U{x,y) = U{x — y) and that D{x,z) = D{z) does not depend on x. Then 
(5.7) is the weak form of, for all x G t>0, 

(5.23) dtf,t{x) = -k D]{x) - y{x)it{x) - a{x)Ctix)[f,tkU]{x), 

where, for example, [j^tkDfix) = JgdCt{dyhiy)D{x - y). 

5.2. Convergence to a superproeess. In this section we show the rela¬ 
tionship between the original BPDL model (rigorously written in Definition 
2.5) and the superprocess version of the Bolker-Pacala model introduced by 
Etheridge [6] . More precisely, we show that accelerating the rates of produc¬ 
tion and natural death by a factor of n makes the BPDL processes converge 
to a continuous random measure-valued process which generalizes the one 
studied in [6]. 


Assumption C2. 
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1. The space T" = The initial conditions Xq converge in law, for the weak 

topology on to a (random) measure Xq G Furthermore, 

sup„-E((X((,l)3) < +00. 

2. There exist some continuous positive functions a{x),a{x),'y{x) and f3{x) 

on respectively bounded by and P, a nonnegative symmetric 

continuous function U{x,y) on x bounded by U, such that 


(5.24) 


ln{x) 

= n-i{x) + (5{x), 

Tn{x) 

= n-f{x). 


a{x) 

Unix) 


n 

Un{x,y) 

II 

Dn{x,z) 

/ n \ 

V27riT(a:) / 


n|zp \ 
2cr(a:)y 


Note that Dn{x,z) is the density of a Gaussian vector with mean 0 and 
variance With these coefficients and when n tends to infinity, we 

have more and more seed production and natural death, and less and less 
competition. Each seed falls more and more close to its mother. 


Theorem 5.6. Admit Assumption C2 and consider the sequence of pro¬ 
cesses X” defined in Notation 5.1. Then for all T > 0, the sequence (X”) 
converges in law, in B([0,T],Mi7’(M'^)), to the unique {in law) superprocess 
X G (^([O,T],Mj?(M'^)), defined by the conditions 

(5.25) sup E[(Xt, 1)^] < oo 

te[o,T] 


and, for any f G (M'^), 

Mt ={XtJ) - {Xo,f)-\ [ ds [ Xs{dx)a{x)-f{x)Af{x) 

Jo JRrf 


fi{x) - a{x) / Xs{dy)U {x, y) 


(5.26) 

- [ ds [ Xs{dx)f{x) 

Jo Jm.'>- 

is a continuous martingale with quadratic variation 

(5.27) {M'^)t = 2 [ ds f Xs{dx)'y{x)f‘^{x) 

Jo JRd 


Proof. We break the proof into several steps. 


Step 1. Let us first prove the uniqueness of the solution of the mar¬ 
tingale problem defined by (5.25)-(5.27); that is, the uniqueness of a prob¬ 
ability measure P on C([0,T],Mi7’(M'^)) under which the canonical process 







22 


N. FOURNIER AND S. MELEARD 


X satisfies (5.25)-(5.27). This result is well known for the super-Brownian 
process (dehned by a similar martingale problem, but with a = f3 = 0 and 
(7 = 7 = 1). As noted in [6], we can use the version of Dawson’s Girsanov 
transform obtained in [5], Theorem 2.3, to deduce the uniqueness in our 
situation, provided the condition 



P{x) 


a{x) / Xs{dy)U{x,y) 


< +00 


is satished. This is easily obtained from (5.25) since the coefficients are 
bounded. 


Step 2. Next we obtain some moment estimates. First we check that 
for all T < oo, 

(5.28) sup sup 1)^] < oo. 

n te[o,r] 

To this end, we use Lemma 5.2(i) with 4){y) = (p, 1)^. [To be completely 
rigorous, first use (/>(p) = (p, 1)^ A A and then make A tend to infinity.] We 
obtain, using Assumption C2, that for all t > 0, all n, 

E [(xr,i)3] 


= E[{X^,lf]+[ ds 
Jo 


E 


Xg{dx)[n^'y{x) +nf3{x)] 


(5.29) 




n 


-(A^1)3 


+ / ds E 


X^{dx)<ri^'y{x)+na{x) / X^{dy)U{x,y) 


1 3 


{X7A)- 


n 


-{x^Af 


Neglecting the nonpositive competition term, we get 

E[{XAl?] 

<E[(X-,1)3] 

XAdx)n‘^j{x) 


+ / ds E 


(5.30) 


n3 


(^r,i) + 


n 


+ 


+ / ds E 


XAdx)nP{x) 


iKA)- 


113 

n 

ll3 


n 


-2{x^Af 
-{x^Af 
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However, for all x > 0, all e G (0,1], (x + e)^ — x^< 6 e(l + x^) and |(x + e)^ + 
(x — e)^ — 2x^1 = 6 e^x. We finally obtain 


(5.31) 


E[{X^, 1)3] < E[{XS, 1)3] + 67 f ds E[{X^, 1)' 

2o 


+ 6/3/ dsE[{X'^,l) + {X^,lf]. 


Part 1 of Assumption C 2 and the Gronwall lemma allow us to conclude that 
(5.28) holds. 

Next, we have to check that 


(5.32) 


supA'f sup (A”,l) j < 00 . 
^ Vte[o,T] / 


Applying Lemma 5.2(iv) with / = 1 and Assumption C2, we obtain 


(5.33) 


{X^,1) = {XS,1)+ f ds [ X^idx) 
Jo JRd. 


fj{x)-a{x) f X2{dy)U{x,y) + . 

JRd J 


Hence 


(5.34) sup {X^, 1) < {XS, 1)+P ds {X2, 1) + sup |M, 


n,l ] 


sG[0,t] 


sg[0,t] 


Thanks to the Doob inequality, part 1 of Assumption C2 and the Gronwall 
lemma, there exists a constant Ct that is not dependent on n such that 

(5.35) e{ sup (A", 1 )) < Ct{l + 

VsG[0,t] / 


Using (5.6) now and Assumption G2, we obtain, for some other constant Ct 
not dependent on n. 


(5.36) £;[(M-- 3 )i ]<(27 + ^) f dsE[{X^,l)]+aU f dsE[{X^,l)^]<Ct 

Jo Jo 

thanks to (5.28). This concludes the proof of (5.32). 


Step 3. We first endow with the vague topology. The exten¬ 

sion to the weak topology is handled in Step 5. We prove the tightness 
of the sequence of laws (£(A"'))„ in V{S}{[0,oo), MpiM.^))) by following 
the same approach as in Theorem 5.3. First, we deduce from Step 2 that 
sup„ i3[sup^g[Q 2 ^] |(A”,/)|] < 00 for any bounded /. We thus have to prove 
that for any / G C'^(M'^), the sequence (Af, /) satisfies the Aldous-Rebolledo 
criterion. Let us consider a couple {S, S') of stopping times satisfying a.s. 
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0<S< S'<S + 6 <T. Using Lemma 5.2, Assumption C2 and the fact 
that I/jjd dz Dn{x, z)f{x + z) — f{x)\ < fj|| A/||oo/2n, we deduce the exis¬ 
tence of a constant C independent of n such that the finite variation part 
of - {XgJ) is bounded by 


(5.37) 


rS+(5 


ds / X^idxMfl 


rS+S 


+ 


ds [ X'^{dx)n^{x) 


dzDn{x,z)f{x + z)- f{x) 


rS+S 


+ / ds X2idx)aU\\f\\^ / X2{dy) 
Js JRd JRd 

rS+S 

<C dsi{X^,l) + {X^,lf). 

Js 


We can also show that, for some constant C, 


(5.38) E[{M'''^)s+s - < CE 


rS+S 


ds{{x^,i)+{x2,iy^ 


Using the moment estimate (5.28), we hnally obtain that the laws of 
and the laws of the drift parts of (A”, /) are tight and then, by Rebolledo’s 
criterion, the laws of (A”,/) are tight. 


Step 4. Let us identify the limit. Let us set Q” = C{X^), denote by 
Q a limiting value of the tight sequence Q” and denote by A = {Xt)t>o a 
process with law Q. Exactly as in the proof of Theorem 5.3, we can show 
that A belongs a.s. to C'([0,T],Mi;’(M'^)). We have to show that A satis¬ 
fies conditions (5.25)-(5.27). First note that (5.25) is straightforward from 
(5.28). Then, we show that for any function / in C'^(M‘^), the process m/ 
dehned by (5.26) is a martingale (the extension to every function in is not 
hard). We consider 0<si<---<Sfc<s<t and some continuous bounded 
maps (/>i,..., on Our aim is to prove that, if the function T from 

D([0,T],Mi7’(M'^)) into M is defined by 

T(z/) = (j)i{iysi) ■ 

(5.39) X |(z/t,/) - (p^,/) du{uu,jaAf/2) 

- [ du [ Uu{dx)f{x) 

Js JR<i 


P{x)- Vu{dy)a{x)U{x,y) 


then 

(5.40) 


E(T(A)) = 0. 
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We know from Lemma 5.2 that using Assumption C2, 

(5.41) 0 = E[MK) • • • = mixn] - An, 

where An is defined by 


An — E 


du / X^{dx) 


(5.42) 


x| 7 (x)n [ dzDnix,z)f{x + z)-f{x)-^^Af{x) 

I iJRd 2n 

+ /3(x) [ dzDnix,z)f(x + z)-f(x) 


First, an easy computation using Assumption C2, that / is and (5.28) 
shows that 

(5.43) |A„| < ^ /* duE[{X::, 1)] ^ 0 

n Js 

as n grows to infinity. Next, it is clear from Assumption C2, the fact that 
/ is cl and that Q only charges the space of continuous processes that the 
map ^ is Q-a.s. continuous. Furthermore, 

(5.44) \^{u)\<c{l +{vs,l) + du{vu,l)‘^'^ 

and we easily deduce from (5.28) that the sequence (|4'(A"')|)„ is uniformly 
integrable. Hence, 

(5.45) \\^E{\^{X'-)\) = Eq{\^{X)\). 

Associating (5.41), (5.43) and (5.45) allows us to conclude that (5.40) holds 
and thus is a martingale. 

We hnally have to show that the bracket of is given by (5.27). To this 
end, we first check that 

Ni = {Xt,ff - {Xo,ff - f ds [ Xs{dx)2j{x)f\x) 

Jo JR'i 


- ds2{XsJ) f Xsidx)fix) 

Jo Jr'^ 

X (5{x)-a{x) [ Xs{dy)U{x,y) 


- / ds2{Xs,f) 


Xs{dx)\a{x)'y{x)Af{x) 


(5.46) 
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is a martingale. This can be done exactly as for m/, using the fact that, 
thanks to Lemma 5.2(iii) (with q = 2), 

W./)" 

- [ ds [ Xg{dx)'y{x) [ dz f‘^{x + z)Dn{x, z) + f‘^{x) 

Jo Jr‘^ L./r'* 

- f ds 2 {X^J) ! K{dx) 

Jo JKd 

X f 3 {x) f dz f{x + z)Dn{x,z) - a{x)f{x) [ X^{dy)U{x,y) 
, , . jRd- JR'' 

(5.47) 

- / ds 2 {X^J) [ X^{dxMx)n 

Jo JR'' 

[ dz f{x +z)Dn{x,z) - f{x) 

.JR'' J 


- [ ds [ X^{dx)f5{x) [ dz p{x + z)Dn{x,z) 

n Jo JR'' JR'' 

-- f ds [ X^{dx)a{x) [ X'^{dy)U{x,y)f{x) 
n Jo JR'' JR'' 

is a martingale for each n. Next, using the Ito formula in the definition (5.26) 
of m/, we deduce that 

{Xtjf-{Xojf-{Mf)t 


(5.48) 



ds2{XsJ) 


Xs{dx)f{x) 



Xs{dy)U{x,y) 


- [ ds2{Xs,f) [ Xs{dx)\a{x)'-f{x)Xf{x) 

Jo JR'' 

is a martingale. Comparing this formula with (5.46) allows us to conclude 
that (5.27) holds. 


Step 5. The extension to the case where MpCU!^) is endowed with the 
weak topology uses similar arguments as in Step 6 of the proof of Theo¬ 
rem 5.3. □ 


6. About extinction and snrvival. First of all, we recall a result in [ 6 ]. 
Consider the superprocess X obtained in Theorem 5.6, and assume that a, 7 , 
P and a are constant on Suppose also that U(x, y) = h{\x — y\) for some 
nonnegative decreasing function h on M_|_ that satisfies /g°° h{r)r‘^~^ dr < 00 . 
Then if /? is sufficiently small and a is sufficiently large, X does not survive: 
a.s., there exists a t > 0 such that for all s > 0 , = 0 . 
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We can also find a complementary result in [6] which shows nonextinction 
with positive probability for another model—the stepping-stone version of 
the Bolker-Pacala process. Let us now come back to the BPDL process 
defined as the solution of (2.8). The techniques used in [6] are specific to 
continuous processes and cannot be generalized to the BPDL discontinuous 
process. 

Before giving our results, let us point out the following obvious remark. 

Remark 6.1. Assume Assumption A and that E[{vq, 1)] < oo. Consider 
the BPDL process {yt)t>o- Assume also that there exist some constants 
7 o < 1*0 such that for all x G A, /r(x) > fiQ and 7 (x) < 70 . Then {vt)t>o does 
a.s. not survive, that is, s > 0, 1) = 0] = 1. 

The proof of this remark is not hard. In such a case, the process Zt = (t't, 1) 
can be bounded from above by a standard continuous-time binary Galton- 
Watson process Yt with death rate /io and birth rate 70 . Since /tq > 70 ) 
extinction a.s. occurs. 

In this section, we first prove almost sure extinction in a case where the 
state space X is compact. Then we show nonextinction in the case of a 
discrete version of the BPDL process with a specific (and not quite realistic) 
competition kernel U. 

6.1. Extinction in the compact case. We check a result which essentially 
says that if the state space X is compact, then the population does almost 
surely not survive. Let us make the following assumption: 

Assumption E. 

(i) The maps a{x) and /r(x) + a{x)U{x,x) are bounded below. 

(ii) There exists a nondecreasing function ip: M_|_ 1 —> M+, satisfying 99 ( 0 ) = 

0, such that limaj^oo ^(x) = 00, such that the map xp{x) is convex on [0, 00) 
and such that, for all n £ Ai, 

(6.1) (u(8)i/,17) > (u, l)(/?((u,l)). 

Remark 6.2. Assumption E(ii) holds if X is compact in and if there 
exist e > 0 and 5 > 0 such that U{x,y) > els^\x-y\<5}- 

Theorem 6.3. Admit Assumptions A and E, i/q G At and E{{uq, 1)) < 00. 
Consider the corresponding unique BPDL process {nt)t>o obtained in Theo¬ 
rem 3.1. Then there is almost surely extinction, that is, P( 31 > 0, (uj, 1) = 
0 ) = 1 . 
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Proof of Remark 6.2. First of all, we cover with a family 
of disjoint cubes of with side Note that L is clearly finite and 

that for each I and each x^y & Ci, \x — y\ <6. Recall the following conse¬ 
quence of the Cauchy-Schwarz inequality, which says that for all L > 1 and 
all {Q;i,...,ai} in R, Hence for all n > 1 and all 

Xi,...,Xn&X, 


n L 


Y. U{Xi,Xj) > 


( 6 . 2 ) 


hi=i 


*j=i 


=-E 


i,j=l 1=1 


1=1 Li=l 


1 2 




Xi 


1 

-^L 


L n 


n 2 




L/=l i=l 


= £—n . 
Jj 


We immediately deduce that for any u € A4, since v is atomic, (zz ®v^U)'> 
e-^(z/, 1)^. Hence Assumption E(ii) holds with ^p{n) = £^n. □ 


Proof of Theorem 6.3. We break the proof into several steps. 


Step 6 . We first of all prove that 


(6.3) A = supE((z't, 1)) <-|-oo. 

t>o 


To this end, we set f{t) = ^((z't,!)) and use Proposition 3.4 with (/)(p) = 
(p, 1) to obtain 


(6.4) f{t) = f{0)+f dsE - y) - [ [ Us{dx)i's{dy)a{x)U{x,y) 

Jo I Jx Jx 


Hence / is differentiable. If we set = ||7 — y\\oo and ao = a{x), we 

deduce that for any t > 0, 


(6.5) f'{t) < 6f{t) - aoE{{ut 0 i^t, U)). 


Using Assumption E and then the Jensen inequality, we obtain that 

(6.6) f'{t)<bf{t)-aof{t)(p{f{t)). 

Let now xq be the greatest solution of bxo = aoXo<f(xo) [recall that ip(x) 
is nondecreasing and goes to infinity with x, and that (/9(0) = 0]. Then we 
deduce from (6.6) that for any f > 0, f(t) < /(O) V xq- This concludes the 
first step. 


Step 7. We now check that a.s. 


(6.7) 


liminf(z't, 1) £ {0, oo}. 
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Since {vt, 1) is N-valued, it suffices to check that for any M G N*, 


P 


liminf(i/t, 1) = M 

>co 


= 0 , 


but this is clear: If liminfi^oo(z^t) 1) = then {vt, 1) reaches the state M 
infinitely often, but reaches the state M — 1 only a finite number of times. 
This is (a.s.) impossible because each time reaches the state M, the 

probability that its next state is M — 1 is bounded below by 


( 6 . 8 ) 


_T/gp_ 

My + Mp, + aUJVP 


> 0 , 


where gp = + a{x)U (x,x)] > 0. 


Step 8. Since {vt, 1) is N-valued and 0 is an absorbing state, we imme¬ 
diately deduce from (6.7) that a.s. lim 4 ^oo(Pt, 1) exists and 

(6.9) lim (pt, 1) G {0, oo}. 

t—^OO 


Step 9. By Fatou’s lemma and Step 1, 


( 6 . 10 ) 


E 


lim {ut, 1) 

t^OO 


= E 


lim inf (pt, 1) 

t—*oo 


< lim inf EUi't, 1)1 < A. 

t-l-OC 


Hence 1) < oo a.s. and we deduce from (6.9) that lim^^oo Wt,l) = 

0 a.s. This concludes the proof. □ 


6.2. Survival in a simplified case. Next, we show that in some cases, 
the BPDL process survives with positive probability. We are not able to 
handle a proof in a general case, because the problem seems very difficult. 
It actually looks much more difficult than the problem of survival for the 
contact process, which has been studied by many mathematicians (see [11]). 
The only result we are able to prove is deduced from a comparison with the 
contact process. 


Assumption S. 


(i) 

(ii) 

(hi) 

(iv) 

( 6 . 11 ) 


The state space X = 

The competition kernel U is pointwise, that is, U{x,y) = l{ 3 ,=y}. 


The dispersion measure Zl(x, dz) = D{dz) = (1/2'^) XlnGZ'^,|n|=i 
7 , pi and a are positive constants that satisfy 


72 


—d 


>2. 


pL + a 


Note that X = was not covered by our construction. The adaptation 
is, however, immediate. 
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Proposition 6.4. Admit Assumption S, assume that vq £ AA, {uq, 1) > 1 
a.s. and assume that EK^q, 1)] < oo. Consider the corresponding BPDL pro¬ 
cess This process survives with positive probability. That means that 

P(inft>o(pt,l) > 1 ) > 0. 

We do not handle a completely rigorous proof. To do so we would have 
to build a rigorous coupling between the contact process and the BPDL 
process. 

Proof or Proposition 6.4. We split the proof into two steps. 


Step 10. Let us first recall definitions and results about the contact 
process (see [11], Chapter VI). First, denote by Mp the set of nonnegative 
finite measures rj on such that for all x £ Z*^, r]{{x}) £ {0, 1}. The contact 
process with parameters > 0 and Am > 0 is a Markov process {T]t)t>o, 
taking its values in Mp and with generator K, defined for all bounded and 
measurable maps (p from into M and all 7j £ Mp{TA) by 


K(p{r]) = XdJ^^r]{dx) ^ l{,,({x+«})=o}['^(?? + 4+«) - </>(??)] 


( 6 . 12 ) 


,\u\=l 


+ Am v{dx)l{n{{x})=i} [4>iv - dx) - 4>iv)]- 


Consider an (possibly random) initial state rjo in Mp satisfying {rjQ,l) > 1 
a.s. Then it is known (see [11], Chapter VI) that the contact process {r]t)t>o 
with parameters \d > 0, Am > 0 and initial state rjQ exists, is unique (in law) 
and that under the condition A^ > 2Am, survives with positive probability. 


Step 11. Consider now the BPDL process (pt)t>o, which takes its values 
in the integer-valued measures on Z'^. Denote pt = J2xe'Z<^ '^{vt(,{x})>i}dx- Note 
that f)t is always dominated by vt. Then {fjt)t>o is a process with values in 
Alp and we can observe that is a sort of contact process with time- 

and space-dependent, random parameters Xd{t, x, w) = 72“'^[1 V Pt({ic})] and 
Xm{t,x,Lo) = lj^j({^}.)<i(/.i-|-a). Under Assumption S, Xd{t,x,ui) is uniformly 
bounded from below by A,^ = 72 “'^, while Xm,{t,x,uj) is uniformly bounded 
from above by Am = fj, + a. Hence, the process {fjt)t>o is bounded below 
by a contact process with parameters Xd Am- Since (6.11) ensures that 
2Am < A^, the conclusion follows from Step 1. □ 


Note that the previously described method may not apply to the continuous- 
state BPDL process, since we really need the interaction to be strictly lo¬ 
cal. In fact, the only case we could treat by such a method is the case 
where the competition kernel is completely loeal and cannot propagate; 
for example, A = and U{x,y) < '^Cp{x)iCp{y), where, for p £ Z'^, 

Cp= [pi,pi + l] X ••• X \pd,Pd + T\- 
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7. On equilibria. An interesting question is that of the existence of non¬ 
trivial equilibria for the BPDL process. Since this question seems very com¬ 
plicated, we first try to give some results about the deterministic equation 

(5.7) . Then we show that there exists a nontrivial equilibrium for the BPDL 
process that is related to the carrying capacity under a detailed balance 
condition which is unfortunately very restrictive. We finally present some 
simulations. We suppose Assumption B in the whole section. 

7.1. Equilibrium of the deterministic equation. We first of all point out 
a trivial remark. 

Remark 7.1. Suppose Assumption B and that 7 < /r, and consider a 
nonnegative finite measure ^0 on Consider the corresponding unique 
solution (^t)t>o C C'([0, 00 ), Mi7’(R'^)) of (5.7). Then tends to 0 as t grows 
to infinity in the sense that 1) < 

This remark follows from a straightforward application of (5.7) with / = 1 
and of the Gronwall lemma. We next generalize the existence of solutions to 

(5.7) to the case of possibly nonintegrable initial conditions. 

Proposition 7.2. Admit Assumption B. Consider a nonnegative bounded 
measurable function ^0 on 

1. There exists a unique function such that: 

(i) for all t>0 and all x G R*^, ^t(x) > 0; 

(ii) for all T < 00 , supig[o,T],a;eR‘* < 0 °; 

(iii) for all t>0 and all x G R*^, 

(7.1) ^t{x) = ^o{x) + f ds['y{^s^D){x)-fi^s{x)-a^s{x){^s*U){x)], 

Jo 

where, for example, (.^t *D)(x) = /jjd dyDix - y)it{y)- 

2. For all x G R*^, the map 1 1 —> f,t{x) is of class on [0,oo), and for all 
T < 00 , |(9t^i(x)| is bounded on [0,r] x R'^. 

3. If furthermore J^d f.o{x) dx < 00 , then for all T < 00 , 

sup / dxf,t{x) < 00 
te[o,r]2Rrf 

and the finite measure-valued function {f^t{x) dx)t>o is the unique solution 
to (5.7). 

Since this proposition is quite unsurprising, we only sketch the proof. 

Proof of Proposition 7.2. First note that point 2 is an immedi¬ 
ate consequence of (7.1) and of the fact that f is bounded, obtained in (i) 
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and (ii). Point 3 is also easily deduced from point 1. To check the unique¬ 
ness part of point 1, it suffices to consider two solutions (?t(a;))i>o 
{Ct{x))t>o xgR'* both bounded by some constant At on [0, T] x 

A straightforward computation shows that, for = supg<j |.^s(x) — 
|s(x)|, for t<T, 


(7.2) 4>{t) < + iJ. + 2aAT) f ds4>{s). 

Jo 

[Recall that since fj^d U{x) dx = 1, sup^j^gd(.^s * U){x) < sup^-g^d Cs(a;).] The 
Gronwall lemma allows us to conclude that ^ 

The existence part follows from an implicit Picard iteration. Define (x) = 
^o(x) and construct by induction a sequence of functions (■Cf)t>o such that 
for each x G 1 1 —> ^t'(x) is of class on M'*' and satisfies, for n > 1, 


(7.3) 


£r+‘(j:)={oW 

+ r * |7(e*0)W-(<{?+'W-a{U'W(C*U(a^)l- 

Jo 


We can, moreover, check at each step that is well dehned, nonnegative and 
bounded on [0, T] x for each n and each T. A straightforward computation 
shows that for all t > 0, sup„ sup^^gj^d (x) < sup 3 ,g]jjd ^o(a^)e'^*) and next that 
for any T, there exists a constant Bt such that for all t <T, 


(7.4) 


sup |{,”+■(!)-{?WI 

XGK'^ 


< Bt 



ds 


sup -c(^)i + sup ic(^) -C“'(^)i 

-xGR'* a;GR‘* 


Thanks to the Gronwall lemma, we deduce that for all T, all t < T and all 
n, 


(7.5) sup|Cr+'(^) 

xGR'* 


cr(a;)l <sBrexp(TRr) f ds sup \C{x) - C Ha:)|. 
Jo xeR'^ 


The Picard lemma allows us to conclude that for all T, 
(7.6) -?r(a;)l <oo- 

„>14G[0,T],xGR‘^ 


Hence, there exists a function ('Ct(x))^>Q ^-gj^d such that for any T, sup^gjQ ^.gjgd (x) — 
^P(x)| tends to 0. We easily check that this function satisfies points (i)-(iii). 

□ 


We may now dehne the equilibria. 




DESCRIPTION OF A REGULATED POPULATION 


33 


Definition 7.3. Admit Assumption B. For a nonnegative bounded con¬ 
tinuous function / on define the function Ff on by 


(7.7) Ff{x) 

Then (7.1) can be rewritten as 


l[f*D]{x) 

^ + a[f-kU]{x)' 


(7.8) ^t{x) = (o{x)+ [ ds{f^i + a[^s*U]{x)){F^s{x)-^s{x)). 

Jo 

This leads us to dehne the equilibria in the following sense. A continuous 
bounded nonnegative function c on is said to be a reasonable equilibrium 
of (7.1) if for all x G 


(7.9) 


c{x) = Fc{x). 


This definition is slightly restrictive, but we may note that if D and U 
are continuous, then any solution to (7.9) such that 

limsup[c*D](x)/[cTt/](x) < oo 

|a:|^oo 

will be continuous and bounded. 


Remark 7.4. Assume Assumption B, that 7 > /.t and that a > 0. Then 
the constant function co(x) = (7 — /r)/a is a reasonable equilibrium of (7.1). 
The constant function c(x) = 0 is also, of course, a reasonable equilibrium 
of (7.1). 

Note that the quantity (7 — /i)/a appears in [2] and is called the earrying 
eapacity, which can be understood as a sort of maximum number of plants 
per unit of volume. We use the following estimate. 

Lemma 7.5. Assume Assumption B, that 7 > /i and that a > 0. Define 
the signed funetion R on by R{x) = D{x) + J^[d{x) — U{x)). Then, for 

all bounded functions f and all x G M'^, 

(7.10) Ff{x) - Fco{x) = rf.Try A if - co)*R]ix). 

p -|- a[f -k U\[x) 

This result is immediately proved by using simply the expression of F. 
We now state an assumption which ensures that R{x) dx is a probability 
measure and hence that T is a contraction around cq in the space of bounded 
functions. 

Assumption C. 7 > /.t and for all x G ^D{x) > (7 — p)U{x). This 
implies that R{x) dx is a probability measure on M'^. 
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Let us now describe a situation for which the constant function cq is the 
unique nontrivial reasonable equilibrium. 

Proposition 7.6. Assume Assumptions B and C, that 7 > 2 '^/r and 
that a > 0. Suppose also that D{x) = D{\x\), where the map D is nonin¬ 
creasing on [0,oo). {This hypothesis is physically reasonable] see [2].) Then 
any nontrivial reasonable equilibrium c of (7.1) identically equals cq. 

Proof. Let c thus be a nontrivial reasonable equilibrium for (7.1). 

Step 1 . Since c is nontrivial, there exists xq such that c { xq ) > 0. Since 
c is continuous, we deduce that c is bounded below on a neighborhood of xq. 
Then (7.9) and the fact that D charges any neighborhood of 0 (since it is 
nonincreasing) ensure that c never vanishes. 


Step 2. We now show that there exists a constant Eq > 0 such that 
for all X G M'^, c{x) > Eq- To this end, we first consider e > 0 such that 
7(1/2'^ — e) > /i and then consider a > 0 such that /[q „]d D{x) dx > 1/2^^ — e, 
which is possible since D is radial. Consider now any point x = (xi,..., x^) G 
and the box B = [xi, xi ->r a] x ■ ■ ■ x [xd, + a]. Denote m = uiix^B c(x), 
which is positive since c is continuous and never vanishes. Our aim is to 
show that m > g{m), where the function g is defined on [ 0 ,oo) by 


(7.11) 


9{u) = f 

f{u) = 


2 ^- 


— £ 


'yu 


// + 07/(7 -/r)u 


This concludes the proof of Step 2 since we can check that g\0) = (1/2'^ — 
£)'y/n > 1 so that m > Eo > 0 where eq is the smallest positive solution to 
u = g{u). 

We thus check that m > g{m). Let y ^ B. Using (7.9) and Assump¬ 
tion C, we deduce that c{y) > f{[c-kD]{y)). However, / is nondecreasing, so 
that c{y) > f{m dzD{y — z)). Using the symmetry and the nonincreas¬ 
ing properties of D, we easily deduce that since y G B, /g dzD{y — z) > 
/[o aY dzD{z) > 1/2^^ — E. Thus for all y G B, c{y) > /(m(l/2'^ — e)) = g{m), 
which ends Step 2. 


Step 3. Using (7.10), Step 2 and Assumption C, we obtain 
sup |c(x) — Col = sup |Fc(x) — Pco(x)| 

(7.12) ^ 

< - sup |c(x) — Co|. 

T + iGRd 

This implies that sup^-gRd |c(x) — co| = 0. □ 
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Although the above uniqueness result seems quite promising, we are at 
the moment not able to prove that under the conditions of the previous 
proposition, any solution {^t)t>o to (7.1) starting from a nontrivial initial 
condition converges to cq in some sense. We can, however, obtain two partial 
results. 

Assumption DBC. a>0, 7>0, ^ = 0 and D = U. 

This assumption is a detailed balance condition. Indeed, under this condi¬ 
tion, the equilibrium co(x) = '^la ensures that for any couple of points x and 
y, the rate of appearance of plants at x due to seed production at y equals 
the rate of disappearance of plants at x because of competition of plants 
at y. In other words, 'yD{x — y)co(y) = aco{x)co{y)U{x — y). Unfortunately, 
this condition is very restrictive. 

Proposition 7.7. Take Assumptions B and DBC. Let be a positive 
bounded and measurable function on Consider the associated unique 
solution of (7.1) starting from obtained in Proposition 7.2. Then 

ft tends to cq = ^la as t grows to infinity in the sense that for all x and all 
t, 

(7.13) [^t(x) - cof < [Co(3:) - co]^exp(-2a[(.fo A cq) -k D]{x)t). 

We furthermore see in the proof below that the behavior of ft is quite 
simple: If fo{x) < cq, then ft{x) increases to cq, while if foix) > cq, then 
ft{x) decreases to cq. 

Proof of Proposition 7.7. Since in this case, dtffix) = —aftkD{x) x 
{ft{x) — Co), we easily show that for all f > 0 and all a; G 

(7.14) dtlffix) - co]^ = -2a[ft{x) - col'^ift* D]{x). 

Since f is nonnegative, we deduce that [ft{x) — cq]^ is nonincreasing in t for 
each x. Since furthermore ft{x) is continuous in t for each x, we deduce that 
for any t,x, ftix) > ^0(2 :) A cq. Hence 

(7.15) dtiftix) -co]^ < -2a[ft{x) - co]^[(Co A cq) *D](a:), 
from which the conclusion follows. □ 

We now treat quite a general case of coefficients a, 7, y, U and D, but 
we consider an initial condition which is only a small perturbation of cq. 

Proposition 7.8. Admit Assumptions B and C, that a > 0 and that 
U is bounded below by a positive continuous function h on 'Mfi. Consider 
a nonnegative bounded measurable function fo on such that /Rd['Co(T) — 
cq]^ dx < 00. Consider the associated unique solution {ft)t>o of (7.1) starting 
from ^0 obtained in Proposition 7.2. Then ft tends to cq as t grows to infinity 
in the sense that there exists a > 0 such that for all t, 
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(7.16) f [6(a^)<e / [Coix] - cof dx. 

JWd J^d 

Proof. We break the proof into three steps. 

Step 1. A straightforward computation using part 2 of Proposition 7.2, 
(7.8) and (7.10) shows that for all f > 0 and all x G 

dt[Ct{x) - cq]"^ 

= 2[?t(a^) - co]dt^t{x) 

= 2[6(a:) - Co] [fJ- + * U){x)] [F^t{x) - 6(a;)] 

= 2[^t{x) -co][ix + a(^f*[/)(x)][F^t(x) -Fco(x)] 

(7.17) +2[^t(a:) - co][ir + a(^t * 17)(a:)][co -^(a;)] 

= 2n[^t{x) -co][(6 -co)*7?](a:) 

- 2[6(a:) - Co]^[// + a(6*77)(x)] 

= -2a[6(a;) - co]^(6*t7)(x) 

- 2//[6(a;) - co][(6(a;) - Co) - {(6 - co)-kR}{x)]. 
Integrating this differential inequality against time, we obtain 

[^i(x) - co]^ 

(7.18) = [^o(a;)-co]^-2 / dsa[(s{x) - cof[^s*U]{x) 

Jo 

-2 [ dsn[^s{x) -co]{[6(a:) -Co] - [{is - cq) -k R]{x)} ds. 

Jo 

Thanks to Assumption C, i? is a probability measure. We furthermore know 
that it, and thus it * U, is bounded on [0,T] x for each T. Thus an 
application of the Cauchy-Schwarz and Young inequalities yields 

(7.19) [ dx [it{x) - co][{itix) - Co)-kR{x)] < f dx [itix) - cof ■ 

jRd jTg^d 

We easily deduce that for all T > 0, supper] /jr^ dx [^t(x) — co]^ < oo. Hence 
(7.18) may be integrated on x G and we get that for all t > 0, 

(7.20) dt [ dx[it{x) - cof <-2a [ dx[it{x) - co]‘^[itkU]{x). 

jRd jRd- 

Step 2. We now wish to bound [ifkU]{x) from below. First, we deduce 
from (7.20) that J^d dx [^t(x) — co]^ is nonincreasing in time. Hence there 
exists a constant b <oo such that for all t > 0, 

dx 1{^^(3 .)<co/2} < b. 


(7.21) 
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However, since U{x) > /i(x), for some positive continuous function h there 
exists a constant a > 0 such that 

(7.22) inf f dzU{z)>ba. 

dx<bJM‘^/A 

Indeed, choose any compact subset K of whose Lebesgue measure equals 
2b and set a = inf^eK h{x). Note that for all A G such that dx < b, 

we also have dx > b, so that 


(7.23) [ dzU{z)> [ dzh{z)>ba. 

JM.'i/A Jk/A 

Finally using (7.22) with A = At^x = {u ^ — y) > co/2}, of which the 

Lebesgue measure is smaller than b thanks to (7.21), we obtain for all x G 
and all t > 0, 


(7.24) 


[^t-kU]{x)= [ dyit{x - y)U{y) 
jR'i 


>co r 

“ 2 Jm,. 


dyU{y) > 


bacQ 


Step 3. Gathering (7.20) and (7.24), we finally obtain 
(7.25) dt [ dx [^i(x) - cof < -bacoa [ dx [^t{x) - cq]^ 

JR'' JR-' 

from which the conclusion follows. □ 


7.2. Equilibrium of the BPDL process. We now to show that it might 
be possible to hnd an equilibrium for the BPDL processes. This is a first 
step to study the long time behavior of the BPDL process {i't)t>o dehned in 
Dehnition 2.5 conditioned on nonextinction. We unfortunately are able to 
treat only the case where the detailed balance condition holds. Of course, 
such an equilibrium will be infinite. We can, however, state the following 
rigorous result. 

We hrst of all denote by M the set of nonnegative (possibly inhnite) 
integer-valued measures on We also denote by A the set of functions (f> 

from A4 into M of the form </>(p) = for some bounded measurable 

function F on M and some function / with compact support on M'^. 

Proposition 7.9. Admit Assumptions B and DBG (see Section 7.1) 
and that U{0) = 0. Consider a Poisson measure vr on with intensity 
measure codx, where CQ = ^la. Then ir is a stationary BPDL process in the 
sense that for all 4> G A, Lfij:) a.s. exists, belongs to and ElLfij:)] = 0, 
where L is defined in (2.3). 
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Note that allowing Assumption DBG and that U{0) = 0 implies that there 
is no natural death. We remark also that this result is somewhat surprising, 
since it suggests that at equilibrium, the plant locations are independent. 
Let us hnally mention that a similar result without Assumption DBG would 
be much more interesting. However, the stationary process tt does not seem 
to be Poisson in such a case. The proof relies on the following lemma, known 
as Slivnyak’s formula in [13] and also can be obtained from Palm measure 
considerations (see [8], Ghapter 10). 


Lemma 7.10. Let u be a Poisson measure on with intensity m{dx). 
Denote by {xi}i>i the points of n, that is, v = Then for all measur¬ 

able functions h from x TW into M such that fj^d m(dx)E[lh(x, v + (Jj,)]] < 
oo, 


(7.26) 


E 


i>l 


h{xi,v) 


= / m{dx)E[h{x,v + 6x)]. 


Proof of Proposition 7.9. Let cj) belong to A. The fact that L(j){7r) 
a.s. exists and belongs to for cf ^ A can be easily checked using the 
explicit expression of L and standard results about Poisson measures. We 
thus prove only that ii^[L(/)(7r)] = 0. Denote by {xi}i>i the points of vr, that 
is, TT = ^xi - Hence, we obtain, using Assumption DBG, 


E[L^{7r)] 



dz D(z){(/>(7r + 6xi+z) - <(>(7r)} 


(7.27) 


- 4i) - </>(vr)} - ■ 


j>i 


+ aE 

.i>l 

=: 7A1 + aA 2 . 

We hrst use Lemma 7.10 with the function hi{x,u) = fj^d dzD(z}{f)(n -h 
dx+z) - 


Ai = E 


(7.28) 


'^hl{Xi,TT) 


i>l 


= codxE / dz D{z){(l){Tr6x + dx+z) - + dx)} ■ 

Jm'i 

Next, with h 2 {x, u) = {(^(p — dx) — </>(i^)} i'{dy)D{x — y), we obtain 

A2 = e( ^h2{xi,7r)\ 


(7.29) 


. i>l 


= f dx cqE {(j){ 7 r) - (j){TT -\-dx)} [ (tt + dx){dy)D{x - y) 
L JR'^ 
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Since -D(O) = U{0) = 0, we obtain, setting h^{ii, v) = D{x 

4)}, 


(7.30) 


Ao = 


dxcoE 



Using Lemma 7.10 again, we obtain 


y){(j){u) - (j){u + 


(7.31) 


A2= dxco dycoE[D{x-y){4>{7r + Sy)-(l){T^ + Sx + Sy)}] 

JM.d. J^d 

= cl[ dx [ dzE[D{z){(j){7r + 6x)-(i)iTr + 6x+z + dx)}], 

JRd- 


where we have used in the last equality the substitution {y,x) i—> (x,x + z). 
Since acg = yco, we deduce that ^Ai = —aA 2 , which ends the proof. □ 


7.3. Simulations. The previous results suggest that the BPDL process, 
conditioned on nonextinction, should converge as time tends to infinity to 
a random measure Voo, quite well distributed (not far from the Lebesgue 
measure), with (7 — /r)/a plants per unit of volume on average. We present 
simulations of this situation. 

We assume that T = M and that 7 = 5, /r = 1 and a = 1. We consider 
the case where U{x,y) = '^{\x-y\<i/2} and D{z) = gl{| 2 |< 3 }- Then we com¬ 
pare the BPDL process {vt)t>o with the stationary solution co{dx) = [(7 — 
y,)/a] dx of (7.1). 

On Figure 1, we assume that t'o = <^o- The boxes represent the empirical 
density of the BPDL process at times t = 3 [Figure 1(a)] and then t = 25 
[Figure 1(b)], obtained with one simulation, while the dotted line is the 
density of cq [i.e., (7 — /r)/a]. We check that after some time, the BPDL 
process is quite well approximated by cq. 

Figure 2 represents the evolution in time of r't([—5,5]) (full line), start¬ 
ing either from uq = 60 [Figure 2(a)] or from uq = 60(5o [Figure 2(b)], and 
compares it with co([—5,5]) = 10(7 — fi)fa (dotted line). 


(a) 


(b) 


jn 


-50 ^0 -30 -20 -10 0 10 20 30 40 50 


-40 -30 -20 -10 0 10 20 30 40 50 


Fig. 1. (a) t = 3; (b) i = 25. 
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(a) (b) 




Fig. 2. (a) 1^0 = So; (b) t'o = 605o. 


Finally, we measure the power of competition. To this end, we compare 
the evolution in time of the rate of interaction of all particles on particles 
located in a ball. We assume that 1^0 = Sq. Figure 3(a) represents, in full line, 
the evolution in time of J^h't{dx) J^i't{dy)l\x\< 5 U{x,y) obtained by one sim¬ 
ulation. The constant value (dotted line) is /j^ co{dx) co{dy)l\x\< 5 U{x, y) = 

10 * [(7 — /i)/a]^. Figure 3(b) shows the same quantities replacing 5 by 50. 

In conclusion, we can say that, on one hand, cq seems to be a good de¬ 
terministic approximation of the BPDL process after a long time. On the 
other hand, there are clearly stochastic fluctuations around the deterministic 
approximation that could be interesting to study. 

Acknowledgments. The authors thank Regis Ferriere and Bernard Roynette 
for numerous fruitful discussions. 
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